Magnetospheric Scattering and Emission in Millisecond Pulsars 
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ABSTRACT 

We model the formation of magnetospheric components of miUisecond pulsar light curves, 
deriving an approximate model for the curved space, 'swept-back' dipole field and following 
photon emission, propagation, and scattering. Magnetospheric pulse components are strongly 
affected by rapid rotation and Schwarzschild effects. 

Subject headings: stars:neutron — pulsars:general — magnetic fields 



1. Introduction and Light Curve Model 

With the launch of CXO and XMM-Newton, 
there has been a dramatic increase in sensitivity 
for high resolution X-ray studies of compact ob- 
jects. Rapidly rotating neutron stars are targets 
of particular interest, since measurements of the 
stellar surface can probe the neutron star evolu- 
tion and the equation of state at high density. Sev- 
eral relativistic effects may also be visible, in some 
cases. A number of millisecond pulsars (MSP) 
have been observed in the X-rays, both as iso- 
lated rotation powered objects (Becker & Triimper 
1999) and as accretion powered LMXBs (e.g. Wij- 
nands & van der Klis 1998). In Braje, Romani, & 
Rauch (2000), we explored the distortions of the 
surface emission imparted by rapid rotation, in- 
cluding the important effects of Doppler boosting, 
aberration, and gravitational focusing and time 
delays. More subtle effects induced by frame drag- 
ging were also considered. For rotation-powered 
pulsars, the surface radiation must traverse the 
magnetosphere, where resonant scattering can in- 
troduce significant distortions in the pulse profile 
(Rajagopal & Romani 1997). Further, several X- 
ray MSP have narrow, non-thermal pulse compo- 
nents, suggesting a direct origin in the magneto- 
sphere. In this paper, we extend the treatment 
of rapid rotation effects to include scattering and 
emission in a surrounding dipole magnetosphere. 

Our model assumes a spherical star, which we 
take to be sufficiently centrally condensed that 



a Schwarzschild (or Kerr) metric is an adequate 
model of the external spacctime (the 'Roche ap- 
proximation'). Because we wish to extend the 
modeling to small spin periods P*, the effects of 
'sweep-back' on the magnetic field, even for r a few 
times the stellar mass M, can be substantial. In 
addition, the rapid rotation develops large electro- 
motive forces. We follow the common assump- 
tion that this EMF causes pair production such 
that the closed zone of the magnetosphere is filled 
with a charge-separated pair plasma, whose charge 
distribution cancels the rotational EMF. For the 
magnetic field, we assume a point dipole located at 
the stellar center and generalize the computation 
of the field structure for the non-rotating dipole 
in the Schwarzschild spacetime (e.g., Bardeen & 
Press 1973; Prasanna & Gupta 1997; Muslimov & 
Harding 1997). Our result passes smoothly to the 
Schwarzschild magnetic field B as the stellar an- 
gular velocity fi* vanishes and to the flat space ro- 
tating dipole solution as r — > oo. We assume that 
the plasma co-rotates (is stationary) in the closed 
zone, which is defined by the 'last closed' field lines 
traced from tangent approach to the light cylinder 
at tlc = cP* /27r to the stellar surface. 

Our light curve modeling follows the proce- 
dure described in Braje et al. (2000). This Monte 
Carlo code starts with photons randomly drawn 
from the surface emission zone, with initial direc- 
tions drawn from a model limb-darkened distri- 
bution. We aberrate these surface photons and 
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propagate to infinity through curved space, us- 
ing the Schwarzschild or Kerr metric, as appro- 
priate. As the photon passes through the mag- 
netosphere, we monitor for local cyclotron reso- 
nance. At the resonance position photons are re- 
emitted with the proper boosted scattering angu- 
lar distribution. We also follow photons arising 
directly from the magnetosphere — from acceler- 
ation gaps or other non- thermal sources. After 
including the gravitational and time-of-flight de- 
lays, the photons are assigned to energy and ro- 
tational phase bins to produce maps of the radia- 
tion on the sky. Slices through these maps at the 
Earth's viewing angle provide pulsar light curves 
and phase-resolved spectra. 

2. Magnetic Field Structure 

We have developed an approximate, retarded, 
dipolar magnetic field expression that links the ex- 
act non-rotating Schwarzschild result at small r 
with the flat space 'swept-back' field structure at 
large r. We expect the magnetic field lines to have 
approximately the same shape as a static dipole 
near the star, but to curve back near the light 
cylinder, as Figure 1 displays. 

To derive our magnetic field expression, we 
write Maxwell's Stress Tensor as derivatives of the 
vector potential in covariant form as 



(1) 



We use the Schwarzschild metric in geometrized 
units {G = c = 1): 



ds^ = rfdt" - ^ - r^d^^ _ ^2 gij^2 Q^^i 



(2) 



where rj = ^J\ — 2M/r. If we solve Maxwell's 
Equations in vacuum 



(3) 



for a dipole moment /j, = we arrive at the so- 
lution in the Schwarzschild spacetime 

/Lt sin^ 6 



A^{r,e) = fir)'- 



(4) 



where we have written the expression as the flat- 
space result weighted by the multiplicative func- 
tion: 

__log,^ + _ l + _ (5) 



(e.g., Wasscrman & Shapiro 1983). For the more 
general rotating case, we multiply the flat space 
time dependent vector potential by the same func- 
tion /(r): 
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H(t)x r U{t)x r 
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where we now have a time dependent dipole mo- 
ment, fi{t) = sinacosl^*(t— r), — sinasinf2^<(^- 
r) , cos a) , tilted with respect to the rotational (z) 
axis by an angle a. We use {et,e^,e6),e<^} to de- 
note the coordinate basis vectors and {e£,ef,eg,e^} 
for the orthonormal basis vectors. By taking the 
dot product of the vector potential with the coor- 
dinate basis vectors, we obtain the vector potential 
components {At, Ar, Aq, A^). Using 



Be = 



Bd, = Fef = 
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(c.f. Muslimov & Harding 1997) and the definition 

of Maxwell's Stress Tensor (eqn. 1), we derive an 
expression for the field: 



Br = 



Ba = 



Bi = 
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(8)' 

with X = 2M/r, x = ^ogif + x{l + x/2), and 

(j) = (f) — Vlf.{t — r) . 

The above expressions recover the familiar re- 
sults in the appropriate limits: for M ^ we ob- 
tain the fiat space retarded fields and for ^ 
we get the Schwarzschild expression for a dipole 
tilted with an angle a with respect to the z 
axis. Checking Maxwell's Equations (3), we find 
that two are exactly satisfied {v = 1,4). The 
other two are nonzero with leading terms of or- 
der 0[{Q.^r/cf{GM/r^)]. 

In Figure 1, we display the field line differences 
induced at small r in curved space coordinates. 
We have drawn field lines from the stellar surface 
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Fig. 1. — Left Panel - Retarded dipolar magnetic field viewed down the rotational axis. The light cylinder 
is depicted as a circle on top of the field lines. Right Panel - Curved space magnetic field perturbations 
viewed in the /i-fi* plane. Flat space (dashed lines) and curved space (solid lines) — the 'last closed' field 
lines are drawn bolder. The points show the closed zone resonant scattering locus for Schwarzschild (circles) 
and flat space (stars) fields, P* = 1.5ms, and the same narrow range of photon energy (2.8eV-3.0eV). B is 
set by the same magnetic moment /Lt(r = oo). 



in the i7«-/Li plane. Note that the field lines are 
pulled closer to the star by curved space distor- 
tions. The open field lines are nearly radial and 
thus display smaller curved space distortions. 

3. Magnetospheric Scattering 

Surface photons {e.g. from a thermal polar cap) 
must propagate through the magnetosphere sur- 
rounding the pulsar. Photons having an energy in 
the correct range pass through a local cyclotron 
resonance in the magnetosphere. For each pho- 
ton energy, this defines a re-emission shell in the 
closed zone (Figure 1) in which strong scattering 
off the re-directs fiux along the local magnetic 
field line direction. For millisecond pulsars with 
B ~ 10^ — 10^ gauss, closed zone resonance typi- 
cally occurs for photons of energy 1 — lOeV. With 
typical surface temperatures of a heated polar cap 
at lO^K, this radiation falls in the Rayleigh-Jeans 
part of the blackbody spectrum. 

In our Monte Carlo code, each photon trajec- 
tory represents a range of emitted photon ener- 
gies. For each photon energy Ej, we find the lo- 
cal cyclotron resonance position along the curved 



space path as the photon redshifts to lower en- 
ergy. For rapidly rotating stars, we must include 
Doppler shifts by boosting the redshifted photon 
energy into the local plasma rest frame, check- 
ing for resonance, re-emitting in this locally co- 
rotating frame, and then boosting back to the lab 
frame where wc continue to integrate the photon 
trajectory. We do not scatter at open zone reso- 
nance positions, since these are assumed to have 
no stationary plasma. 

In principle, photons can undergo multiple scat- 
terings in the magnetosphere, for sufficiently large 
optical depth t. Following Rajagopal & Romani 
(1997), we estimate r assuming a stationary, flat 
space dipolar magnetic field. The cyclotron energy 
as a function of distance from the star is 

= 11.6Bi*2(l + 3cos2 6'B)^/^(i?*/r)3keV (9) 

where BI2 is the equatorial surface dipole field 
strength in units of 10^^ gauss, 6b is the polar an- 
gle between the dipole moment and the direction 
of the local magnetic field, and i?* is the neutron 
star radius. The resonant cyclotron cross section 
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where e_ contains the dependence on the pho- 
ton direction and polarization and ap is the fine 
structure constant (e.g., Meszaros 1992). We as- 
sume a co-rotation (Goldreich- Julian) charge den- 
sity riGj = 7 X 10^^Bf2/-P*(ms) cm"^ (Goldreich 
& Julian 1969). For rapid magnetospheric pair 
production the total density may, of course, 
be higher. Integrating the resonant cross section 
through the co-rotation charge density we get the 
optical depth for scattering 



n(r)c7resdr 



2.3 

P* (ms) 



B 



.1/3 /^3eVy/'B 



B lOkih 



where we have taken = 1/3. 

For concreteness, we adopt fiducial parameters 
for a millisecond pulsar (P* = 3ms, P* = 10km, 
and P* = 10" gauss); these give a characteristic 
optical depth of r = 7.7 x IQ-^B^/B for E^ = 
3eV. For polar cap emission, the average optical 
depth is typically somewhat smaller as significant 
flux travels through the scattering-free open zone. 
Thus, for co-rotation charge densities, we neglect 
multiple scatterings. 

In Figure 2, we plot light curves for thermal 
emission from a single polar cap, including closed 
zone scattering. The asyinnicitry in the scattered 
flux for the simple Schwarzschild propagations is 
due to the field 'sweep-back'. Doppler boosts, 
aberration, and time delay effects are quite impor- 
tant, producing asymmetry in the direct flux and 
further phase shifting the scattered components. 
The 'inter-pulse' peaks in the scattered flux are 
due to radiation scattered back past the star and 
gravitationally focused to produce peaked pulse 
components. The strong energy dependence of 
the location of the scattering screen can been seen 
in the shifts of this pulse component. Note that 
the scattering amplitude decreases with Ej as is 
apparent from equation (11). The energy depen- 
dence of the phase shifts is larger for rapid rotators 
(small P*), as the high altitude scattering acquires 
additional aberration and Doppler boosting. For- 
mally, these scattered radiation components pro- 
vide a precision probe of the near star gravita- 
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Fig. 2. — Normalized MSP light curves for a sin- 
gle thermal cap and viewing angle ( = 115°, in- 
cluding closed zone resonant scattering. Left — 
P* = 1.5ms. Right - P* = 3.0ms. Upper panels - 
simple Schwarzschild metric propagations. Lower 
panels - propagations include Doppler shifts and 
time delays. Curves are plotted for E^ = 1.6eV 
(top) to 16eV (bottom). 



tional field, with field geometry and other systcm- 
atics calibrated through the energy dependence 
of the component phase. In principle, precision 
tests of strong gravity are possible. Wc have run 
the code with Kerr metric propagations to check 
the detectability of these effects. At P* = 1.5ms 
we find frame dragging shifts the positions of the 
scattered photon peaks (Figure 2) by At/) ~ 0.01. 
We conclude that unless neutron stars are very 
compact at short periods, or unless closed zone 
plasmas substantially exceed the co-rotation den- 
sity, frame dragging effects will be quite difficult to 
measure in scattered photons from surface thermal 
emission. 

4. Outer Gap Emission 

As a straightforward extension of these mod- 
els, we consider the effects of spacetime curvature 
on the light curves from outer gap emission, up- 
dating the work of Romani & Yadigaroglu (1995) 
to include curved space propagations. The outer 
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gap region follows the edge of the closed zone ex- 
tending from the 'null charge surface' (where ncj 
changes sign, i.e. Bz = 0) to tlc- In outer gap 
models this region remains charge-starved, large 
potential drops develop, particles accelerate, and 
non-thermal radiation is produced tangent to the 
local B (e.g., Cheng, Ho, & Ruderman 1986; Ro- 
mani 1996). For MSP, this emission zone comes 
very close to the neutron star surface where space- 
time curvature becomes important. 

To find the null charge boundary, we make a 
bicubic spline of the closed zone surface and step 
out along the grid until _Bz = 0. We emit tangen- 
tially from the 'last closed' field line surface with 
an emissivity F cx r~^/^. As described in Romani 
& Yadigaroglu (1995) aberration, boosts, and time 
delays are all essential for any outer gap light curve 
computation. Accordingly, here we check the ad- 
ditional distortions induced by curved space ef- 
fects, as might be important for millisecond pul- 
sars. 
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Fig. 3. — Light curves for outer gap emission. 
P* = 1.5ms, a = 65°, C = 115°. The magnetic 
axis (polar cap) is at phase 0.5. Schwarzschild 
light curves (solid line) differ from flat space curves 
(dashed line) in both pulse phase and width. 

In Figure 3, we show sample light curves for 
outer gap emission. The phase of the magnetic 
axis may be inferred from the structure of the ra- 
dio polarization data and so the relative phasing of 
the high energy (outer gap) pulse is significant. Of 
course even classical pulsar radio emission, if aris- 
ing from more than a few i?*, will show similar 
phase shifts, so comparison of light curve phases 



must account for the altitude of the radio emis- 
sion zone. The effect of neglecting Schwarzschild 
propagations is to both compress the pulse and 
shift it to later phase. For our fiducial parameters 
the first pulse component arises at high altitude. 
The second pulse arises from small r where the 
cap field lines (and photon directions) are nearly 
radial, allowing only small gravitational bending. 
The curved space differences are thus dominated 
by the phase delays induced by the retardation in 
the strong near surface field. The curved space 
pulse expansion (relative to flat space) can sim- 
ilarly be attributed to gravitational time delay, 
since the first peak arises at larger r where these 
eff'ects are small. 

5. Conclusions and Observational Prospects 

We have improved the modeling of light curves 
of rapidly rotating neutron stars (MSP) by ex- 
tending the description of the surrounding vac- 
uum dipole magnetosphere to include curved space 
and rapid rotation effects up to corrections of or- 
der 0\{Vl^r/cf{GM/rc^)]. These effects produce 
substantial changes in the pulse components in- 
troduced by magnetospheric resonant scattering 
(for thermal surface emission) or direct emission 
from the rotating magnetosphere (high altitude 
gap radiation). We have developed a Monte Carlo 
code to compute light curves illustrating these ef- 
fects including the energy dependence of the scat- 
tered photon pulse shapes. Higher order eflfects 
from frame dragging have also been computed, but 
these are likely too subtle to be discernible in mag- 
netospheric pulse components unless the S/N is 
very high. Scattering perturbations, though small, 
will be greatly enhanced if the magnetosphere can 
support a plasma with r > tgj- In this paper, 
we have only considered the minimum possible 
scattering perturbation. With new detector tech- 
nologies (e.g. Romani et al. 1999), we should see 
a marked increase in the quality of optical light 
curves. With enhanced sensitivity and spectral 
resolution (and possibly higher scattering r), light 
curve perturbations due to rapid rotation effects 
may be promoted to important probes of pulsar 
physics. 
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grant NAG5-3263. Roger W. Romani is a Cottrell 



5 



Scholar of Research Corporation. 
REFERENCES 

Bardeen, J.M. & Press, W.H. 1973, Journal of 
Mathematical Physics, 14, 7 

Becker, W. & Triimper, J. 1999, A&A, 341, 803 

Braje, T.M., Romani, R.W., & Ranch, K.P. 2000, 
ApJ, 531, 447 

Cheng, K.S., Ho, C, & Ruderman, M. 1986, ApJ, 

300, 500 

Coldrcich, P. & Julian, W.H. 1969, ApJ, 157, 869 

Meszaros, P. 1992, High Energy Radiation from 
Magnetized Neutron Stars (Chicago:University 
of Chicago Press) 

Muslimov, A. & Harding, A. K. 1997, ApJ, 485, 
735 

Prasanna, A. R. & Gupta, A. 1997, Nuovo Ci- 
mento B, 112B, 1089 

Rajagopal, M. & Romani, R.W. 1997, ApJ, 491, 
296 

Romani, R.W. 1996, ApJ, 470, 469 

Romani, R.W., Miller, A.J., Cabrera, B., 
Figueroa-Feliciano, E., & Nam, S.W. 1999, 
ApJ, 521, 153 

Romani, R.W. & Yadigaroglu, I.-A. 1995, ApJ, 
438, 314 

Wasserman, I., & Shapiro, S. L. 1983, ApJ, 265, 

1036 

Wijnands, R., & van der Klis, M. 1998, Nature, 
394, 344 



This 2-column preprint was prepared with the AAS lAT^^X 
macros v5.0. 



6 



